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Abstract. Stellar winds emitted by rotating massive stars may show two main components: firstly bipolar lobes 
with low density and fast wind, produced by the higher T c fr and gravity at the poles ("g e ff-effect"); secondly, an 
equatorial disc with a slow dense wind, produced by the stronger opacities at the equator ("^-effect"). To see the 
possible role of this anisotropic wind on the shape of LBV nebulae, we calculate the distribution of the ejected 
matter in 2 simplified cases: 1) A brief shell ejection. We find that prolate and peanut-shaped hollow nebulae 
naturally form due to the g c ff-effect in rotating stars. 2) A constant wind for a long time. This produces prolate 
filled nebulae, with a possible strong disc when a bi-stability limit is crossed in the equatorial region. Thus, many 
features of the rj Carinae and LBV nebulae are accounted for by the anisotropic ejection from rotating stars. 
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1. Introduction 

The HST picture of the rj Carinae nebula (J. Morse and K. 
Davidson, STScI PRC96-23a) shows two big symmetrical 
lobes and a dis k-sha ped skirt around t he sta r (Davidson 
& Humphreys |l997| Davidson et al. [l997| ). These au- 
thors find that the lobes were created during the "Great 
Eruption" of 1843, during which ~ 1 to 3 M were ejected 
forming the Homunculus Nebula, with a total energy of 
~ 10 49 5 ergs (Humphreys 1999| ). T he kin ematics confirms 
that it is a bipolar outflow (Nota 1999), and the study 
of their limb darkening shows that the lobes are hollow 
(Davidson et al. 1997| ). The mass in the skirt amounts to 



nario of an aspherical fast wind, expanding into a previ- 
ously deposited slow spherical wind has been studied by 



0.1 to 0.2 M Q , it is still unclear whether it is resulting from 
the 1843 or later events. A massive cold torus of about 15 
Mq has been found out in the equatorial plane (Morris 
et al. 1999). Peanut-shaped nebulae and lobes are present 
around sever al other LBV stars (Nota |1999| ). 

Spectral variations with a period of 5.5 yr have led to 
the suggestion that rj Carinae is a binary system (Damineli 
199C; Lamers et al. 1998; Damineli et al. 2000| ). However, 
the periodicity is not as strict as claimed, with shifts of 3 
monthes between two successive cycles, which casts some 
doubt on the binary hypothesis (Smith et al. |2000D 

Numerous models for the rj Carinae nebula have been 



proposed (e.g. see Schulte-Ladbeck 1997; Hillier 1997). 
Several ones invoke collisions of winds emitted at various 
evolutionary stages, such as the interaction of an isotropic 
fast wind with a previous slow equatorially enhanced wind 
(Kwok et al. |1978| Frank et al. |1995p. The opposite sce- 



Frank et al. (1998). This model recovers several features, 
but not the equatorial skirt. Langer (1998), Langer & 
Heger (1998) consider that the LBV outbursts occur when 
the outwards centrifugal and radiation forces cancel grav- 
ity at the equator. Langer et al. (1999) develop a model 
where a slow, dense, equatorial wind is followed by a fast 
and almost spherical wind and this leads to very represen- 
tative simulations of the Homunculus Nebula. The models 
are using the wind compressed disc model of Bjorkman 
and Cassinelli ( 1993 ). 

Owocki et al. ( 1993 ) hav e emp hasized that the account 
of the von Zeipel theorem ( 1924 ) drastically modifies the 
wind distribution around a rotating star. The properties 
of stellar wind s ejected by rotating stars have been studied 
(Maeder [l99§ Maeder & Meynet |2000| ), in particular the 
changes of the mass loss rates by rotation and the latitu- 
dinal distribution of the mass loss. Here, we examine the 
distribution of the anisotropic wind around rotating star 
models corresponding to r\ Carinae. The noticeable result 
is that there are bi-polar components and a disc natu- 
rally arising in the wind. Interactions of the lobes and 
discs emitted at various phases during the back and forth 
motions of an LBV in the HR diagram may further shape 
the nebulae (cf. Langer et al. 1999[ ). 
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2. The anisotropic mass loss 

In rotating stars with high radiation pressure, the surface 
shape is also described by the Roche model. The total 
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Fig. 1. The mass fluxes around a rotating star like 77 
Carinae for a ratio of the angular velocity to the break-up 
angular velocity u> — 0.80, assuming a polar T e a = 30000 
K. 

gravity is the sum of the gravitational, centrifugal and 
radiative accelerations: g to t = 9cS + fifrad = Sgrav + fifrot + 
<7rad- The local flux at colatitude # is proportional to the 
effective gravity g c ff according to the von Zeipel theorem 
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where we ignore a small term (~ 1 %) due to differential 
rotation. The correct Eddington factor r^i?) in a rotating 
star depends on rotation (Maeder & Meynet 2000). 
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where is the local opacity. One has 

fftot =fleff[l-r O (0)] ■ 



(2) 



(3) 



The relation g tot — gives the critical velocity. The ex- 
pression t^ rit = — T) often used for massive stars 
only applies if the star is uniformly bright. If not, Eq. 



(3) has two roots which need to be discussed (Maeder & 
Meynet 2000). The result is that, due to the von Zeipel 
theorem, w cr i t only slowly tends towards zero for high 
Eddington factor. 

The distribution of the mass loss by unit surface at 
colatitude $ over the surface of a rotating star is 
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where Tq is given by Eq. (2) for electron scattering opac- 

ity. We have A = (ka) c " ( — —) 01 . The parameter a is a 
so-called force multiplier, which characterizes the intensi- 
ties and distribution of the strengths of the spectral lines. 
At some T e ff, the ionisation equilibrium of the stellar wind 
is changing rather abruptly and so does the opacity of 
the plasma. Consequently, the values of the force multi- 
pliers undergo rapid transitions, e.g. a = 0.54 for log T ef j 
> 4.342, a = 0.30 for log T cff between 4.342 and 4.00 and 
a = 0.20 for log T cff < 4.00 (Lamers et al. |l995| ; Lamers 
1997; see also Kudritzki & Puis | 20o"by Such transitions of 



the wind properties, which still have some uncertainties, 
are called bi-stability limits. Since T e g is decreasing from 
pole to equator, one has the corresponding variations of 
the force multipliers and thus of the mass loss fluxes from 
pole to equator. 

Eq. (4) shows that two main effects may enhance the 
mass ejection. The first one is the higher gravity at the 
pole which makes it hotter, this is the "<? e ff _ effect" . The 
second is the "re-effect", it is due to the rapid growth of 
the opacity below a bi-stability limit. The result is a lower 
value of a in the equatorial regions, if they become cool 
enough; this favours strong equatorial ejections. 

We consider the case of 77 Carinae, with a mass of 
~100 M©, a luminosity of ~ 10 6 5 Lq, an Eddington 
factor r — 0.823 for electron scattering opacity, and a 
T e fj ~ 30000 K or a few thousands degr ees le ss at mini- 
mum (Humphreys 198S; Davidson et al. 1997 ). Rotation 
is defined by u>, the fraction of the angular velocity at 
break-up, i.e. ui 2 = . Fig. 1 shows the distribution 

of the mass fluxes for u> — 0.80 and a reference polar 
T ff = 30000 K. We notice the prolate distribution with 
a peanut-shape, which results from the g e ff-effect. For 
higher values of u, the peanut-shape is more and more 
pronounced, while for lower w we have a prolate spheroid. 

When a bi-stability limit is crossed at some latitude 
on the rotating star, the M-rates are strongly enhanced 
at lower latitudes by the re-effect. A crossing of the bi- 
stability limit is favoured by fast rotation and a lower 
choice of the reference T c ff . Fig. 2 shows, as an example, 
the mass flux for u> = 0.8 as in Fig. 1, but for T e ff = 
25000 K. We notice the presence of a strong equatorial 
ejection. The wider the latitude range with T e ff below the 
bistability limit, the thicker the disk. 
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Fig. 2. The mass flux around a rotating star with the 
mass and luminosity of r\ Carinae for to = 0.80 and a 
polar T cff = 25000 K. 

3. The shape of the nebulae 

The shape of the nebulae is given by the distribution of 
matter at some time around the ejecting star. For purpose 
of simplicity, we consider here two simplified cases: a) a 
short time ejection, b) a continuous ejection. 



3.1. Short time ejection 

We assume that the geometry r(t, d) at time t of the shell 
ejected at to is determined by the coasting velocity, while 
Eq. (4) determines the amount of matter at a given posi- 
tion in the shell. Thus, one has (cf. Maeder |1999| ) 



r(t,0) = uoo(r,0)(t-to) 



with 
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Fig. 3 shows the distribution of matter for a short time 
ejection for a case with u> = 0.90. If there is no bi-stability 
limit crossed, this shape is the same for all kinds of stars 
with a given ui (cf. Eq. 6). The nebula is hollow and shows 
a peanut shape with an equatorial pinch. For lower val- 
ues of u, the shape is a prolate spheroid; for higher u), 
the equatorial pinch is even more marked. At a given ui, 
the peanut shape is more pronounced for the M distribu- 
tion than for the Voa distribution, since the first one de- 

i 

pends on g e g, while the second one depends on g c 2 ff . Such 
a sha pe we ll corresponds to that of many LBV nebulae 
(Nota |1999D. 



Fig. 3. Spatial distribution of the matter for a short ejec- 
tion around a star with u — 0.90. The vertical line shows 
the direction of the rotation axis. 

3.2. Continuous ejection 

We may also assume the wind to be constant in time, and 
for a long time. On each radial line, the density decreases 
like r~ 2 . Thus, the lines r($) of constant densities are 
given by 
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This model gives centrally filled nebulae. Fig. 4 illustrates 

the iso-density lines in the case corresponding to Fig. 1. 

i 

Since, the iso-density lines behave like g e 4 ff , the variations 
of g c s have smaller effects and the peanut shape is absent. 
We notice a strong equatorial enhancement due to the 
lower a in these regions. In real cases, if a sudden ejec- 
tion is superposed on a more continuous wind, we may 
have some compound of the hollow peanut shell as in Fig. 
3 and of the filled nebula with a disc as in Fg. 4, well 
corresponding to the various possible shapes of the LBV 
nebulae and rj Carinae. Interestingly enough, the group 
of the B[e] stars (Zickgraf 1999) shows a two-component 
stellar wind with a hot, highly ionized, fast wind at the 
poles and a slow dense disk-like wind at the equator. 

Figs. 3 and 4 only give two schematic examples of neb- 
ulae. Certainly, a kind of catalogue is to be made of all 
possible sorts of nebulae arising from Eq. (4). This is a 
problem with many parameters: 
- Stellar rotation given by ui, 
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Fig. 4. Iso-density lines for a continuous ejection around 
a star with u> = 0.80. The density is strongly enhanced in 
the equatorial plane. 



Reference polar T e g of the star, 
Eddington factor T, 

Opacities, bi-stability limits, opacity peaks, etc. 
Particular history of mass ejection M(t). 
Wind-wind interaction (cf. Langer et al. 



1999| ) 



Certainly many interesting situations may occur. For ex- 
ample, a B-star with a polar T e g not so much above the 
bi-stability limit could give rise to a substantial equato- 
rial ejection and a disc formation, even if the star is not 
at break-up ! 

Finally, we may wonder whether LBV stars have fast 
rotation velocities. We see two possibilities, a) If due to 
a small magnetic field rotation keeps solid body (Langer 
1998), for most non-zero initial velocities the star may 
reach the break-up velocity near the end of the MS. b) 
Models with differential rotation also show that very mas- 
sive stars may reach the so-called fir-limit, i.e the break- 
limit function of rotation and of the Eddington factor 
(Maeder & Meynet 2000| ). In particular, after the end of 
the MS phase, the stars do a bluewards hook in the HR 
diagram; also after an outburst, they may reach fast ve- 
locities as a result of the contraction of the outer layers. 



4. Conclusions 

Rotation produces strongly anisotropic mass loss in mas- 
sive stars. Two schematic cases of winds have been con- 
sidered. The occurence of a peanut-shaped nebula with 
bi-polar lobes and an equatorial disc may naturally result 
from the the g e g and K-effects. Without particular wind 
interactions and collisions, the wind distribution already 



possesses several important features of 77 Carinae and of 
other LBV stars. 

The critical question is whether the polar lobes and the 
skirt in rj Carinae result directly from the anisotropic mass 
ejection (cf. Figs. 3 & 4) or whether they result from the 
interaction of the successive anisotropic winds. Whatever 
the answer, the anisotropic ejection with polar lobes and 
an equatorial disc, as shown here, is an essential part of 
the game. 

Acknowledgements. A.M. expresses his best thanks to Dr. 
Norbert Langer and to Dr. Georges Meynet for their very help- 
ful and constructive comments, and to Raoul Behrend for his 
help in numerical modelling. 



References 

Bjorkman, J.E., Cassinelli, J.R 1993, ApJ 409, 429 
Damineli, A. 1996, ApJ 460, L49 

Damineli, A., Kaufer, A., Wolf, B., Stahl, O. Lopes, D.F., de 

Arajo, F.X. 2000, ApJ 528, L101 
Davidson, K., Humphreys, R. M. 1997, ARA&A 35, 1 
Davidson, K., Ebbets, D., Johansson, S., Morse, J. A., Hamann, 

F. W. 1997, AJ 113, 335 
Davidson, K., Gull, T.R., Humphreys, R.M., Ishibashi, K., 

Whitelock, P. et al. 1999, AJ 118, 1777 
Frank, A., Balick, B., Davidson, K. 1995, ApJ 441, L77 
Frank, A., Ryu, D., Davidson, K. 1998, ApJ 500, 291 
Hillier J. 1997, ASP Conf. Ser. 120, p. 287 
Humphreys, R. M. 1992, in Physics of Luminous Blue 

Variables, IAU Coll. 113, Rds. K. Davidson et al., Kluwer, 

p.3 

Humphreys, R. M. 1999, Lecture Notes in Physics 523, 243 
Kudritzki, R.P., Puis, J. 2000, ARAA, 38, 613 
Kwok, S., Purton, C.R., Fitzgerald, P.M. 1978, ApJ 210, L125 
Lamers H. 1997, ASP. Conf. Ser. 120, 83 

Lamers, H., Livio, M., Panagia, N., Walborn, N. 1998, ApJ 
505, L131 

Lamers H., Snow TP., Lindholm D.L. 1995, ApJ 455, 269 
Langer, N. 1998, A&A 329, 551 
Langer, N., Heger, A. 1998, ASP Conf. Ser. 131, 76 
Langer, N., Garcia- Segura, G., Mac Low, M.M. 1999, ApJ 520, 
L49 

Maeder, A. 1992, in "Instabilities in evolved super- and hy- 

pergiants", Ed. C. de Jager and H. Nieuwenhuijzen, North 

Holland, p. 138. 
Maeder, A. 1999, A&A 347, 185 
Maeder, A., Meynet, G. 2000, A&A 361, 159 
Meynet, G., Maeder, A. 2000, A&A 361, 101 
Morris, P.W., Waters, L.B.F.M., Barlow, M.J., Lim T. et al. 

1999, Nature, 402, 502 
Owocki, SR., Gayley, K.G., Cranmer, SR. 1998, ASP Conf. 

Ser. 131, 237 

Nota, A., Clampin, M. 1999, ASP Conf. Ser. 120, p. 303 
Nota, A. 1999, Lecture Notes in Physics, 523, 62 
Schulte-Ladbeck R.E. 1997, ASP Conf. Ser. 120, p. 355 
von Zeipel H. 1924, MNRAS 84, 665 

Smith, N., Morse, J.A., Davidson, K. Humphreys, R.M. 2000, 
AJ 120, 920 

Zickgraf, F. J. 1999, Lecture Notes in Physics, 523, 40 



